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Abstract. Engineers in the manufacturing industries have used accel- 
erated test (AT) experiments for many decades. The purpose of AT 
experiments is to acquire rehabihty information quickly. Test units of 
a material, component, subsystem or entire systems are subjected to 
higher-than-usual levels of one or more accelerating variables such as 
temperature or stress. Then the AT results are used to predict life 
of the units at use conditions. The extrapolation is typically justified 
(correctly or incorrectly) on the basis of physically motivated models 
or a combination of empirical model fitting with a sufficient amount of 
previous experience in testing similar units. The need to extrapolate 
in both time and the accelerating variables generally necessitates the 
use of fully parametric models. Statisticians have made important con- 
tributions in the development of appropriate stochastic models for AT 
data [typically a distribution for the response and regression relation- 
ships between the parameters of this distribution and the accelerating 
variable(s)], statistical methods for AT planning (choice of accelerat- 
ing variable levels and allocation of available test units to those levels) 
and methods of estimation of suitable reliability metrics. This paper 
provides a review of many of the AT models that have been used suc- 
cessfully in this area. 

Key words and phrases: Reliability, regression model, lifetime data, 
degradation data, extrapolation, acceleration factor, Arrhenius rela- 
tionship, Eyring relationship, inverse power relationship, voltage-stress 
acceleration, photodegradation. 
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1. INTRODUCTION 

1.1 Motivation 

Today's manufacturers face strong pressure to de- 
velop new, higher-technology products in record time, 
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while improving productivity, product field reliabil- 
ity and overall quality. This has motivated the de- 
velopment of methods like concurrent engineering 
and encouraged wider use of designed experiments 
for product and process improvement. The require- 
ments for higher reliability have increased the need 
for more up-front testing of materials, components 
and systems. This is in line with the modern quality 
philosophy for producing high-reliability products: 
achieve high reliability by improving the design and 
manufacturing processes; move away from reliance 
on inspection (or screening) to achieve high reliabil- 
ity, as described in Meeker and Hamada (1995) and 
Meeker and Escobar (2004). 

Estimating the failure-time distribution or long- 
term performance of components of high-reliability 
products is particularly difficult. Most modern prod- 
ucts are designed to operate without failure for years. 
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decades or longer. Thus few units will fail or degrade 
appreciably in a test of practical length at normal 
use conditions. For example, the design and con- 
struction of a communications satellite may allow 
only eight months to test components that are ex- 
pected to be in service for 10 or 15 years. For such 
applications, Accelerated Tests (ATs) are used in 
manufacturing industries to assess or demonstrate 
component and subsystem reliability, to certify com- 
ponents, to detect failure modes so that they can 
be corrected, to compare different manufacturers, 
and so forth. ATs have become increasingly impor- 
tant because of rapidly changing technologies, more 
complicated products with more components, higher 
customer expectations for better reliability and the 
need for rapid product development. There are diffi- 
cult practical and statistical issues involved in accel- 
erating the life of a complicated product that can fail 
in different ways. Generally, information from tests 
at high levels of one or more accelerating variables 
(e.g., use rate, temperature, voltage or pressure) is 
extrapolated, through a physically reasonable statis- 
tical model, to obtain estimates of life or long-term 
performance at lower, normal levels of the acceler- 
ating variable(s). 

Statisticians in manufacturing industries are of- 
ten asked to become involved in planning or ana- 
lyzing data from accelerated tests. At first glance, 
the statistics of accelerated testing appears to in- 
volve little more than regression analysis, perhaps 
with a few complicating factors, such as censored 
data. The very nature of ATs, however, always re- 
quires extrapolation in the accelerating variable(s) 
and often requires extrapolation in time. This im- 
plies critical importance of model choice. Relying 
on the common statistical practice of fitting curves 
to data can result in an inadequate model or even 
nonsense results. Statisticians working on AT pro- 
grams need to be aware of general principles of AT 
modeling and current best practices. 

The purpose of this review paper is to outline 
some of the basic ideas behind accelerated testing 
and especially to review currently used AT modeling 
practice and to describe the most commonly used 
AT models. In our concluding remarks we make ex- 
plicit suggestions about the potential contributions 
that statisticians should be making to the develop- 
ment of AT models and methods. We illustrate the 
use of the different models with a series of examples 
from the literature and our own experiences. 



1.2 Quantitative versus Qualitative 
Accelerated Tests 

Within the reliability discipline, the term "acceler- 
ated test" is used to describe two completely differ- 
ent kinds of useful, important tests that have com- 
pletely different purposes. To distinguish between 
these, the terms "quantitative accelerated tests" (Qua- 
nAT) and "qualitative accelerated tests" (QualAT) 
are sometimes used. 

A QuanAT tests units at combinations of higher- 
than-usual levels of certain accelerating variables. 
The purpose of a QuanAT is to obtain information 
about the failure-time distribution or degradation 
distribution at specified "use" levels of these vari- 
ables. Generally failure modes of interest are known 
ahead of time, and there is some knowledge avail- 
able that describes the relationship between the fail- 
ure mechanism and the accelerating variables (either 
based on physical/chemical theory or large amounts 
of previous experience with similar tests) that can 
be used to identify a model that can be used to 
justify the extrapolation. In this paper, we describe 
models for QuanATs. 

A QualAT tests units at higher-than-usual combi- 
nations of variables like temperature cycling and vi- 
bration. Specific names of such tests include HALT 
(for highly accelerated life tests), STRIFE (stress- 
life) and EST (environmental stress testing). The 
purpose of such tests is to identify product weak- 
nesses caused by flaws in the product's design or 
manufacturing process. Nelson (1990, pages 37-39) 
describes such tests as "elephant tests" and outlines 
some important issues related to QualATs. 

When there is a failure in a QualAT it is necessary 
to find and carefully study the failure's root cause 
and assess whether the failure mode could occur in 
actual use or not. Knowledge and physical/chemical 
modeling of the particular failure mode is useful for 
helping to make this assessment. When it is deter- 
mined that a failure could occur in actual use, it 
is necessary to change the product design or man- 
ufacturing process to eliminate that cause of fail- 
ure. Nelson (1990, page 38) describes an example 
where a costly effort was made to remove a high- 
stress-induced failure mode that never would have 
occurred in actual use. 

Because the results of a QualAT are used to make 
changes on the product design or manufacturing pro- 
cess, it is difficult, or at the very least, very risky to 
use the test data to predict what will happen in nor- 
mal use. Thus QualATs are generally thought of as 
being nonstatistical. 
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1.3 Overview 

The rest of this paper is organized as follows. 
Section 2 describes the basic physical and practi- 
cal ideas behind the use of ATs and the charac- 
teristics of various kinds of AT data. Section 3 de- 
scribes the concept of a time-transformation model 
as an accelerated failure-time model, describes some 
commonly used special cases and also presents sev- 
eral nonstandard special cases that are important 
in practice. Section 4 describes acceleration models 
that are used when product use rate is increased 
to get information quickly. Sections 5 and 6 explain 
and illustrate the use of temperature and humid- 
ity, respectively, to accelerate failure mechanisms. 
Section 7 describes some of the special characteris- 
tics of ATs for photodegradation. Section 8 explains 
and illustrates the use of increased voltage (or volt- 
age stress) in ATs with the commonly used inverse 
power relationship. This section also describes how 
a more general relationship, based on the Box-Cox 
transformation, can be used in sensitivity analyses 
that help engineers to make decisions. Section 9 de- 
scribes examples in which combinations of two or 
more accelerating variables are used in an AT. Sec- 
tion 10 discusses some practical concerns and gen- 
eral guidelines for conducting and interpreting ATs. 
Section 11 describes areas of future research in the 
development of accelerated test models and the role 
that statisticians will have in these developments. 

2. BASIC IDEAS OF ACCELERATED TESTING 

2.1 Different Types of Acceleration 

The term "acceleration" has many different mean- 
ings within the field of reliability, but the term gen- 
erally implies making "time" (on whatever scale is 
used to measure device or component life) go more 
quickly, so that reliability information can be ob- 
tained more rapidly. 

2.2 Methods of Acceleration 

There are different methods of accelerating a reli- 
ability test: 

Increase the use rate of the product. This method 
is appropriate for products that are ordinarily not 
in continuous use. For example, the median life of 
a bearing for a certain washing machine agitator is 
12 years, based on an assumed use rate of 8 loads per 
week. If the machine is tested at 112 loads per week 
(16 per day), the median life is reduced to roughly 



10 months, under the assumption that the increased 
use rate does not change the cycles to failure distri- 
bution. Also, because it is not necessary to have all 
units fail in a life test, useful reliability information 
could be obtained in a matter of weeks instead of 
months. 

Increase the intensity of the exposure to radia- 
tion. Various types of radiation can lead to material 
degradation and product failure. For example, or- 
ganic materials (ranging from human skin to mate- 
rials like epoxies and polyvinyl chloride or PVC) will 
degrade when exposed to ultraviolet (UV) radiation. 
Electrical insulation exposed to gamma rays in nu- 
clear power plants will degrade more rapidly than 
similar insulation in similar environments without 
the radiation. Modeling and acceleration of degra- 
dation processes by increasing radiation intensity is 
commonly done in a manner that is similar to accel- 
eration by increasing use rate. 

Increase the aging rate of the product. Increasing 
the level of experimental variables like temperature 
or humidity can accelerate the chemical processes of 
certain failure mechanisms such as chemical degra- 
dation (resulting in eventual weakening and failure) 
of an adhesive mechanical bond or the growth of a 
conducting filament across an insulator (eventually 
causing a short circuit). 

Increase the level of stress {e.g., amplitude in tem- 
perature cycling, voltage, or pressure) under which 
test units operate. A unit will fail when its strength 
drops below applied stress. Thus a unit at a high 
stress will generally fail more rapidly than it would 
have failed at low stress. 

Combinations of these methods of acceleration are 
also employed. Variables like voltage and tempera- 
ture cycling can both increase the rate of an electro- 
chemical reaction (thus accelerating the aging rate) 
and increase stress relative to strength. In such sit- 
uations, when the effect of an accelerating variable 
is complicated, there may not be enough physical 
knowledge to provide an adequate physical model 
for acceleration (and extrapolation) . Empirical mod- 
els may or may not be useful for extrapolation to 
use conditions. 

2.3 Types of Responses 

It is useful to distinguish among ATs on the basis 
of the nature of the response. 
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Accelerated Binary Tests (ABTs). The response 
in an ABT is binary. That is, whether the prod- 
uct has failed or not is the only reliability informa- 
tion obtained from each unit. See Meeker and Hahn 
(1977) for an example and references. 

Accelerated Life Tests (ALTs). The response in 
an ALT is directly related to the lifetime of the 
product. Typically, ALT data are right-censored be- 
cause the test is stopped before all units fail. In 
other cases, the ALT response is interval-censored 
because failures are discovered at particular inspec- 
tion times. See Chapters 2~10 of Nelson (1990) for 
a comprehensive treatment of ALTs. 

Accelerated Repeated Measures Degradation Tests 
(ARMDTs). In an ARMDT, one measures degra- 
dation on a sample of units at different points in 
time. In general, each unit provides several degrada- 
tion measurements. The degradation response could 
be actual chemical or physical degradation or per- 
formance degradation (e.g., drop in power output). 
See Meeker and Escobar (1998, Chapters 13 and 21) 
for examples of ARMDT modeling and analysis. 

Accelerated Destructive Degradation Tests 
(ADDTs). An ADDT is similar to an ARMDT, ex- 
cept that the measurements are destructive, so one 
can obtain only one observation per test unit. See 
Escobar, Meeker, Kugler and Kramer (2003) for a 
discussion of ADDT methodology and a detailed 
case study. 

These different kinds of ATs can be closely re- 
lated because they can involve the same underlying 
physical/chemical mechanisms for failure and mod- 
els for acceleration. They are different, however, in 
that different kinds of statistical models and analy- 
ses are performed because of the differences in the 
kind of response that can be observed. 

Many of the underlying physical model assump- 
tions, concepts and practices are the same for ABTs, 
ALTs, ARMDTs and ADDTs. There are close re- 
lationships among models for ABT, ALT, ARMD 
and ADD data. Because of the different types of 
responses, however, the actual models fitted to the 
data and methods of analysis differ. In some cases, 
analysts use degradation-level data to define failure 
times. For example, turning ARMDT data into ALT 
data generally simplifies analysis but may sacrifice 
useful information. An important characteristic of 
all ATs is the need to extrapolate outside the range 
of available data: tests are done at accelerated con- 
ditions, but estimates are needed at use conditions. 



Such extrapolation requires strong model assump- 
tions. 

3. STATISTICAL MODELS FOR 
ACCELERATION 

This section discusses acceleration models and some 
physical considerations that lead to these models. 
For further information on these models, see Nelson 
(1990, Chapter 2) and Meeker and Escobar (1998, 
Chapter 18). Other useful references include Smith 
(1996), Section 7 of Tobias and Trindade (1995), 
Sections 2 and 9 of Jensen (1995) and Klinger, Nakada 
and Menendez (1990). 

Interpretation of accelerated test data requires mod- 
els that relate accelerating variables like tempera- 
ture, voltage, pressure, size, etc. to time acceler- 
ation. For testing over some range of accelerating 
variables, one can fit a model to the data to describe 
the effect that the variables have on the failure- 
causing processes. The general idea is to test at high 
levels of the accelerating variable(s) to speed up fail- 
ure processes and extrapolate to lower levels of the 
accelerating variable(s). For some situations, a phys- 
ically reasonable statistical model may allow such 
extrapolation. 

Physical acceleration models. For well-understood 
failure mechanisms, one may have a model based on 
physical/ chemical theory that describes the failure- 
causing process over the range of the data and pro- 
vides extrapolation to use conditions. The relation- 
ship between accelerating variables and the actual 
failure mechanism is usually extremely complicated. 
Often, however, one has a simple model that ade- 
quately describes the process. For example, failure 
may result from a complicated chemical process with 
many steps, but there may be one rate-limiting (or 
dominant) step and a good understanding of this 
part of the process may provide a model that is ad- 
equate for extrapolation. 

Empirical acceleration models. When there is lit- 
tle understanding of the chemical or physical pro- 
cesses leading to failure, it may be impossible to 
develop a model based on physical/chemical theory. 
An empirical model may be the only alternative. An 
empirical model may provide an excellent fit to the 
available data, but provide nonsense extrapolations 
(e.g., the quadratic models used in Meeker and Es- 
cobar, 1998, Section 17.5). In some situations there 
may be extensive empirical experience with particu- 
lar combinations of variables and failure mechanisms 
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and this experience may provide the needed justifi- 
cation for extrapolation to use conditions. 

In the rest of this section we will describe the gen- 
eral time-transformation model and some special ac- 
celeration models that have been useful in specific 
applications. 

3.1 General Time- Transformation Functions 

A time-transformation model maps time at one 
level of X, say X{/, to time at another level of x. 
This can be expressed as T(x) = T[r(x[/),x], where 
xjj denotes use conditions. To be a time transfor- 
mation, the function T(t,x) must have the following 
properties: 

• For any x, T(0,x) = 0, as in Figure 1. 

• T(t,x) is nonnegative, that is, T(t,x) > for all 
t and X. 

• For fixed x, T(t,x) is monotone increasing in t. 

• When evaluated at x^, the transformation is the 
identity transformation [i.e., T{t,xij) = t for all t]. 

A quantile of the distribution of r(x) can be de- 
termined as a function of x and the corresponding 
quantile of the distribution of T(x[/). In particu- 
lar, tp(x) = T[tp(x[/),x] for < p < 1. As shown in 
Figure 1, a plot of T{xu) versus r(x) can imply a 
particular class of transformation functions. In par- 
ticular, 

• r(x) entirely below the diagonal line implies ac- 
celeration. 

• T(x) entirely above the diagonal line implies de- 
celeration. 

• T(x) can cross the diagonal, in which case the 
transformation is accelerating over some times and 
decelerating over other times. In this case the c.d.f.'s 
of r(x) and T{xij) cross. See Martin (1982) and 
LuValle, Welsher and Svoboda (1988) for further 
discussion of time-transformation models. 

3.2 Scale-Accelerated Failure-Time 
Models (SAFTs) 

A simple, commonly used model used to char- 
acterize the effect that explanatory variables x = 
(xi, . . . , Xk)' have on lifetime T is the scale-accelerated 
failure-time (SAFT) model. The model is ubiquitous 
in the statistical literature where it is generally re- 
ferred to as the "accelerated failure-time model." It 
is, however, a very special kind of accelerated failure- 
time model. Some of the explanatory variables in x 
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Fig. 1. General failure-time transformation with Xu <x. 

are used for acceleration, but others may be of inter- 
est for other reasons (e.g., for product design opti- 
mization decisions). Under a SAFT model, lifetime 
at X, r(x), is scaled by a deterministic factor that 
might depend on x and unknown fixed parameters. 
More specifically, a model for the random variable 
T(x) is SAFT if T(x) =T(xc/)/^.F(x), where the 
acceleration factor AJ-{x) is a positive function of x 
satisfying AT{x.u) = 1. Lifetime is accelerated (de- 
celerated) when AJ-{'x) > 1 [AJ-{yi) < !]• In terms 
of distribution quantiles, 

ip(xc/) 



(1) 



tp(x) 



Some special cases of these important SAFT models 
are discussed in the following sections. 

Observe that under a SAFT model, the probabil- 
ity that failure at conditions x occurs at or before 
time t can be written as Pr[T(x) <t\ = Pr[r(x[/) < 
AJ-{x.) xt\. It is common practice (but certainly not 
necessary) to assume that lifetime T(x) has a log- 
location-scale distribution, with parameters {fi,cr), 
such as a lognormal distribution in which /x is a func- 
tion of the accelerating variable(s) and a is constant 
(i.e., does not depend on x). In this case, 

log(t) 



F{t;xu)=PT[T{xu)<t] = <^ 



a 



where $ denotes a standard cumulative distribution 
function (e.g., standard normal) and /if/ is the lo- 
cation parameter for the distribution of log[T(x[/)]. 
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Thus, 



F(t;x) =Pr[r(x) <t] 



log(t)-{/ic/-log[^^(x)]} 



Note that T(x) also has a log-location-scale distribu- 
tion with location parameter /j, = fifj — log[^^(x)] 
and a scale parameter a that does not depend on x. 

3.3 The Proportional Hazard Regression Model 

For a continuous cdf F(t;x[/) and \I'(x) > define 
the time transformation 

r(x) = F-Hl - {1 - F[r(xf/);xc/]}i/*('');xc/). 

It can be shown that r(x) and T(x[/) have the pro- 
portional hazard (PH) relationship 



(2) 



/i(t;x) = ^'(x)/i(i;x 



U) 



This time-transformation function is illustrated in 
Figure 1. In this example, the amount of accelera- 
tion (or deceleration), T(x{/)/r(x), depends on the 
position in time and the model is not a SAFT. If 
F(t;xij) has a Weibull distribution with scale pa- 
rameter ijjj and shape parameter /?{/, then T(x) = 
r(x[/)/^jr(x), where AT{yi) = [^'(x)]^/'^^^ . This im- 
plies that this particular PH regression model is also 
a SAFT regression model. It can be shown that the 
Weibull distribution is the only distribution in which 
both (1) and (2) hold. Lawless (1986) illustrates this 
result nicely. 
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Fig. 2. Superalloy fatigue data with fitted log-quadratic 
Weibull regression model with nonconstant a. Censored ob- 
servations are indicated by >. The response, cycles, is shown 
on the X-axis. 



3.4 Another Non-SAFT Example: The 
Nonconstant a Regression Model 

This section describes acceleration models with 
nonconstant a. In some lifetime applications, it is 
useful to consider log-location-scale models in which 
both and a depend on explanatory variables. The 
log-quantile function for this model is 

log[tp(x)]=/x(x) + $-l(p)cT(x). 



Thus 



tp(x) 



+ cl>-i(p)[a(x^)-a(x)]}. 



Because tp{xij) / tp{x) depends on p, this model is 
not a SAFT model. 

Example 1 ( Weibull log-quadratic regression mo- 
del with nonconstant a for the superalloy fatigue 
data). Meeker and Escobar (1998, Section 17.5) 
analyze superalloy fatigue data using a Weibull model 
in which /i = + /jj^'x + pI^^x^ and log(a) = /3f ' + 
P^^x (see Nelson, 1984 and 1990, for a similar anal- 
ysis using a lognormal distribution). Figure 2 shows 
the log-quadratic Weibull regression model with non- 
constant a fit to the superalloy fatigue data. 

Meeker and Escobar (1998) indicate that the ev- 
idence for nonconstant a in the data is not strong. 
But having a decrease with stress or strain is typi- 
cal in fatigue data and this is what the data points 
plotted in Figure 2 show. Thus, it is reasonable to 
use a model with decreasing a in this case, even in 
the absence of "statistical significance," especially 
because assuming a constant sigma could lead to 
anti-conservative estimates of life at lower levels of 
stress. 

4. USE-RATE ACCELERATION 

Increasing the use rate can be an effective method 
of acceleration for some products. Use-rate accelera- 
tion may be appropriate for products such as electri- 
cal motors, relays and switches, paper copiers, print- 
ers, and home appliances such as toasters and wash- 
ing machines. Also it is common practice to increase 
the cycling rate (or frequency) in fatigue testing. 
The manner in which the use rate is increased may 
depend on the product. 
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4.1 Simple Use-Rate Acceleration Models 

There is a basic assumption underlying simple use- 
rate acceleration models. Useful life must be ade- 
quately modeled with cycles of operation as the time 
scale and cycling rate (or frequency) should not af- 
fect the cycles-to-failure distribution. This is reason- 
able if cycling simulates actual use and if the cycling 
frequency is low enough that test units return to 
steady state after each cycle (e.g., cool down). 

In such simple situations, where the cycles-to-failure 
distribution does not depend on the cycling rate, we 
say that reciprocity holds. This implies that the un- 
derlying model for lifetime versus use rate is SAFT 
where ^^(UseRate) = UseRate/UseRatec/ is the fac- 
tor by which the test is accelerated. For example, 
Nelson (1990, page 16) states that failure of rolling 
bearings can be accelerated by running them at three 
or more times the normal speed. Johnston et al. 
(1979) demonstrated that the cycles-to-failure of elec- 
trical insulation was shortened, approximately, by a 
factor of AT{A12) = 412/60 « 6.87 when the applied 
AC voltage in endurance tests was increased from 
60 Hz to 412 Hz. 

ALTs with increased use rate attempt to simulate 
actual use. Thus other environmental factors should 
be controlled to mimic actual use environments. If 
the cycling rate is too high, it can cause reciprocity 
breakdown. For example, it is necessary to have test 
units (such toaster) "cool down" between cy- 
cles of operation. Otherwise, heat buildup can cause 
the cycles-to-failure distribution to depend on the 
cycling rate. 

4.2 Cycles to Failure Depends on Use Rate 

Testing at higher frequencies could shorten test 
times but could also affect the cycles-to-failure dis- 
tribution due to specimen heating or other effects. 
In some complicated situations, wear rate or degra- 
dation rate depends on cycling frequency. Also, a 
product may deteriorate in stand-by as well as dur- 
ing actual use. Reciprocity breakdown is known to 
occur, for example, for certain components in copy- 
ing machines where components tend to last longer 
(in terms of cycles) when printing is done at higher 
rates. Bowling (1993, page 706) describes how in- 
creased cycle rate may affect the crack growth rate 
in per cycle fatigue testing. In such cases, the empiri- 
cal power-rule relationship ^^(UseRate) = 
(UseRate/UseRatCfy)^ is often used, where p can be 
estimated by testing at two or more use rates. 



Example 2 (Increased cycling rate for low-cycle 
fatigue tests). Fatigue life is typically measured in 
cycles to failure. To estimate low-cycle fatigue life 
of metal specimens, testing is done using cycling 
rates typically ranging between 10 Hz and 50 Hz 
(where 1 Hz is one stress cycle per second), depend- 
ing on material type and available test equipment. 
At 50 Hz, accumulation of 10^ cycles would require 
about five hours of testing. Accumulation of 10'' cy- 
cles would require about two days and accumulation 
of 10^ cycles would require about 20 days. Higher 
frequencies are used in the study of high-cycle fa- 
tigue. 

Some fatigue tests are conducted to estimate crack 
growth rates, often as a function of explanatory vari- 
ables like stress and temperature. Such tests gener- 
ally use rectangular compact tension test specimens 
containing a long slot cut normal to the centerline 
with a chevron machined into the end of the notch. 
Because the location of the chevron is a point of 
highest stress, a crack will initiate and grow from 
there. Other fatigue tests measure cycles to failure. 
Such tests use cylindrical dog-bone-shaped speci- 
mens. Again, cracks tend to initiate in the narrow 
part of the dog bone, although sometimes a notch is 
cut into the specimen to initiate the crack. 

Cycling rates in fatigue tests are generally increased 
to a point where the desired response can still be 
measured without distortion. For both kinds of fa- 
tigue tests, the results are used as inputs to engi- 
neering models that predict the life of actual system 
components. The details of such models that are ac- 
tually used in practice are usually proprietary, but 
are typified, for example, by Miner's rule (e.g., page 
494 of Nelson, 1990) which uses results of tests in 
which specimens are tested at constant stress to pre- 
dict life in which system components are exposed to 
varying stresses. Example 15.3 in Meeker and Es- 
cobar (1998) describes, generally, how results of fa- 
tigue tests on specimens are used to predict the re- 
liability of a jet engine turbine disk. 

There is a danger, however, that increased tem- 
perature due to increased cycling rate will affect 
the cycles-to-failure distribution. This is especially 
true if there are effects like creep-fatigue interaction 
(see Bowling, 1993, page 706, for further discus- 
sion). In another example, there was concern that 
changes in cycling rate would affect the distribu- 
tion of lubricant on a rolling bearing surface. In 
particular, if T is life in cycles and T has a log- 
location-scale distribution with parameters (/i,0"), 
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then /X = /?o + log (cycles/unit time) where f3o and 
Pi can be estimated from data at two or more values 
of cycles/unit time. 

5. USING TEMPERATURE TO ACCELERATE 
FAILURE MECHANISMS 

It is sometimes said that high temperature is the 
enemy of reliability. Increasing temperature is one 
of the most commonly used methods to accelerate a 
failure mechanism. 

5.1 Arrhenius Relationship for Reaction Rates 

The Arrhenius relationship is a widely used model 
to describe the effect that temperature has on the 
rate of a simple chemical reaction. This relationship 
can be written as 

(3) ^(*^-P)=^oexp(^^^^^) 

where TZ is the reaction rate, and tempK = 
temp °C + 273.15 is thermodynamic temperature in 
kelvin (K), k is either Boltzmann's constant or the 
universal gas constant and Ea is the activation en- 
ergy. The parameters Ea and 70 are product or ma- 
terial characteristics. In applications involving elec- 
tronic component reliability, Boltzmann's constant 
k = 8.6171 X 10"^ = 1/11605 in units of electron- 
volt per kelvin (eV/K) is commonly used and in this 
case, Ea has units of electronvolt (eV). 

In the case of a simple one-step chemical reac- 
tion, Ea would represent an activation energy that 
quantifies the minimum amount of energy needed to 
allow a certain chemical reaction to occur. In most 
applications involving temperature acceleration of a 
failure mechanism, the situation is much more com- 
plicated. For example, a chemical degradation pro- 
cess may have multiple steps operating in series or 
parallel, with each step having its own rate constant 
and activation energy. Generally, the hope is that 
the behavior of the more complicated process can 
be approximated, over the entire range of tempera- 
ture of interest, by the Arrhenius relationship. This 
hope can be realized, for example, if there is a single 
step in the degradation process that is rate-limiting 
and thus, for all practical purposes, controls the rate 
of the entire reaction. Of course this is a strong as- 
sumption that in most practical applications is im- 
possible to verify completely. In most accelerated 
test applications, it would be more appropriate to 
refer to Ea in (3) as a quasi- activation energy. 



5.2 Arrhenius Relationship 
Time-Acceleration Factor 

The Arrhenius acceleration factor is 

^jr(temp, tempj;, Ea) 

^ ^(temp) 
TZite-mpjj) 

\^ ( 11605 11605 \1 

= exp Ea — - — . 

I Vtemp{/K tempK/ 

When temp > temp^, ^jr(temp, temp^, > 1. 
When tempf; and Ea are understood to be, respec- 
tively, product use temperature and reaction-specific 
quasi activation energy, ^.7^(temp) = ^jF(temp, 
tempjj, Ea) will be used to denote a time-acceleration 
factor. The following example illustrates how one 
can assess approximate acceleration factors for a 
proposed accelerated test. 

Example 3 {Adhesive-bonded power element). 
Meeker and Hahn (1985) describe an adhesive-bonded 
power element that was designed for use at temp = 
50 °C. Suppose that a life test of this element is to 
be conducted at temp = 120 °C. Also suppose that 
experience with this product suggested that Ea can 
vary in the range Ea = 0.4 eV to Ea = 0.6 eV. Fig- 
ure 3 gives the acceleration factors for the chemical 
reaction when testing the power element at 120 °C 
and quasi-activation energies of Ea = (0.4,0.5,0.6) 
eV. The corresponding approximate acceleration fac- 
tors at 120 °C are ^^"(120) = 12.9,24.5, and 46.4, 
respectively. 

The Arrhenius relationship does not apply to all 
temperature acceleration problems and will be ad- 
equate over only a limited temperature range (de- 
pending on the particular application). Yet it is sat- 
isfactorily and widely used in many applications. 
Nelson (1990, page 76) comments that ". . . in certain 
applications (e.g., motor insulation), if the Arrhe- 
nius relationship . . . does not fit the data, the data 
are suspect rather than the relationship." 

5.3 Eyring Relationship Time-Acceleration 
Factor 

The Arrhenius relationship (3) was discovered by 
Svante Arrhenius through empirical observation in 
the late 1800s. Eyring (e.g., Gladstone, Laidler and 
Eyring, 1941, or Eyring, 1980) gives physical theory 
describing the effect that temperature has on a re- 
action rate. Written in terms of a reaction rate, the 
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Fig. 3. Time- acceleration factors as a function of 
temperature for the adhesive-bonded example with 
= (0.4,0.5,0.6) eV. 



Eyring relationship is 

7?.(temp) = 70 X A(tenip) x exp 



-Ea 



k X temp K 



where >l(teiiip) is a function of temperature depend- 
ing on the specifics of the reaction dynamics and 70 
and Ea are constants (Weston and Schwarz, 1972, 
e.g., provides more detail). Applications in the liter- 
ature have typically used ^(temp) = (temp K)™ with 
a fixed value of m ranging between m = (Boc- 
caletti et al., 1989, page 379), m = 0.5 (Klinger, 
1991), to m = 1 (Nelson, 1990, page 100, and Mann, 
Schafer and Singpurwalla, 1974, page 436). 

The Eyring relationship temperature acceleration 
factor is 



^J^Ey (temp, temp^, Ea) 
_ / tempK N"" 



temp^ K / 
AJ^Ari'temp, temp^, ^a) 



X 



where ^^Ar(temp,tempf/,£'a) is the Arrhenius ac- 
celeration factor from (4). For use over practical 
ranges of temperature acceleration, and for practi- 
cal values of m not far from 0, the factor outside the 
exponential has relatively little effect on the acceler- 
ation factor and the additional term is often dropped 
in favor of the simpler Arrhenius relationship. 

Example 4 (Eyring acceleration factor for a met- 
allization failure mode). An accelerated life test 
will be used to study a metallization failure mecha- 
nism for a solid-state electronic device. Experience 



with this type of failure mechanism suggests that the 
quasi-activation energy should be in the neighbor- 
hood of Ea = 1.2 eV. The usual operating junction 
temperature for the device is 90 °C. The Eyring ac- 
celeration factor for testing at 160 °C, using m = 1, 
is 



>lJ^Ey(160,90,1.2) 

/160-F 273.15 



X ^j^Ar(160,90,1.2) 



V 90-^273.15 
= 1.1935 X 491 = 586 

where ^^Ar(160,90, 1.2) = 491 is the Arrhenius ac- 
celeration factor. We see that, for a fixed value of 
Ea, the Eyring relationship predicts, in this case, 
an acceleration that is 19% greater than the Arrhe- 
nius relationship. As explained below, however, this 
figure exaggerates the practical difference between 
these models. 

When fitting models to limited data, the estimate 
of Ea depends strongly on the assumed value for 
m (e.g., or 1). This dependency will compensate 
for and reduce the effect of changing the assumed 
value of m. Only with extremely large amounts of 
data would it be possible to adequately separate 
the effects of m and Ea using data alone. If m can 
be determined accurately on the basis of physical 
considerations, the Eyring relationship could lead 
to better low-stress extrapolations. Numerical evi- 
dence shows that the acceleration factor obtained 
from the Eyring model assuming m known, and es- 
timating Ea from the data, is monotone decreasing 
as a function of m. Then the Eyring model gives 
smaller acceleration factors and smaller extrapola- 
tion to use levels of temperature when m > 0. When 
m < 0, Arrhenius gives a smaller acceleration fac- 
tor and a conservative extrapolation to use levels of 
temperature. 

5.4 Reaction-Rate Acceleration for a Nonlinear 
Degradation Path Model 

Some simple chemical degradation processes (first- 
order kinetics) might be described by the following 
path model: 



(5) 



Vi^t; temp) 



^ 'Doo X {1 - exp[-TZu X AT{temp) x t]} 



where TZu is the reaction rate at use temperature 
tempj;, TZjj X AJ-{temp) is the rate reaction at a 
general temperature temp, and for temp > tempj^. 
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^^(temp) > 1. Figure 4 shows this function for fixed 
TZjj, Ea and Voo, but at different temperatures. Note 
from (5) that when Pqo > 0, 'D{t) is increasing and 
failure occurs when V^t) > T>{. For the example in 
Figure 4, however, Voo < 0, 'D(t) is decreasing, and 
failure occurs when T>{t) <T>f = —0.5. In either case, 
equating ViT; temp) to T>f and solving for failure 
time gives 



At any other temperature. 



(6) 



where 



T(temp) 



r(temp^) 
AJ^ (temp) 



is failure time at use conditions. Faster degradation 
shortens time to any particular definition of failure 
(e.g., crossing Df or some other specified level) by 
a scale factor that depends on temperature. Thus 
changing temperature is similar to changing the units 
of time. Consequently, the time-to-failure distribu- 
tions at tempjj and temp are related by 



(7) 



Pr[r(temp^) < t] 

= Pr[r(temp) < t/^j^(temp)]. 



Equations (6) and (7) are forms of the scale-accele- 
rated failure-time (SAFT) model introduced in Sec- 
tion 3.2. 

With a SAFT model, for example, if T(tempf;) 
(time at use or some other baseline temperature) 
has a log-location-scale distribution with parameters 
Hu and o", then 



Pr[r < t; temp^] = ^> 



log(t) - 



a 



Pr[r<t;temp] = $ 



log(t) - ^1 



a 



where 



/i = /i(x) = nu - log[^J^(temp)] = Po + Pix, 

X = 11605/(tempK), xu = 11605/(temp^ K), pi = 
Ea and (3o = f^u — [3ixu- LuValle, Welsher and Svo- 
boda (1988) and Klinger (1992) describe more gen- 
eral physical/chemical degradation model charac- 
teristics needed to assure that the SAFT property 
holds. 

Example 5 {Time acceleration for Device- A). 
Hooper and Amster (1990) analyze the temperature- 
accelerated life test data on a particular kind of 
electronic device that is identified here as Device-A. 
The data are given in Meeker and Escobar (1998, 
page 637). The purpose of the experiment was to 
determine if Device-A would meet its failure-rate 
objective through 10,000 hours and 30,000 hours 
at its nominal operating ambient temperature of 
10 °C. Figure 5 shows the censored life data and 
the Arrhenius-lognormal ML fit of the distribution 
quantiles versus temperature, describing the rela- 
tionship between life and temperature. There were 
failure out of 30 units tested at 10 °C, 10 out of 100 
at 40 °C, 9 out of 20 at 60 °C, and 14 out of 15 at 
80 °C. The censored observations are denoted in Fig- 
ure 5 by A. The life-temperature relationship plots 
as a family of straight lines because temperature is 
plotted on an Arrhenius axis and life is plotted on a 
log axis. The densities are normal densities because 
the lognormal life distributions are plotted on a log 
axis. 



^6- 




Fig. 4. Nonlinear degradation paths at different 
tures with a SAFT relationship. 




tempera- FiG. 5. Arrhenius-lognormal model fitted to the Device-A 
data. Censored observations are indicated by A. 
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5.5 Examples where the Arrhenius Model is 
not Appropriate 

As described in Section 5.1, strictly speaking, the 
Arrhenius relationship will describe the rate of a 
chemical reaction only under special circumstances. 
It is easy to construct examples where the Arrhenius 
model does not hold. For example, if there is more 
than one competing chemical reaction and those 
chemical reactions have different activation energies, 
the Arrhenius model will not describe the rate of the 
overall chemical reaction. 

Example 6 {Acceleration of parallel chemical re- 
actions). Consider the chemical degradation path 
model having two separate reactions contributing to 
failure and described by 

V^t; temp) 

= T^ioo X {1 — exp[—TZiu X ^:Fi(temp) x t]} 
+ V200 

X {1 — exp[— 7^21/ X AT2{te-m.-p) x t]}. 

Here IZiu and IZ2U are the use-condition rates of the 
two parallel reactions contributing to failure. Sup- 
pose that the Arrhenius relationship can be used 
to describe temperature dependence for these rates, 
providing acceleration functions .AJ^i(tenip) and 
^.7^2 (temp). Then, unless A!Fi{temp) = ^jF2(temp) 
for all temp, this degradation model does not lead 
to a SAFT model. Intuitively, this is because tem- 
perature affects the two degradation processes dif- 
ferently, inducing a nonlinearity into the accelera- 
tion function relating times at two different temper- 
atures. 

To obtain useful extrapolation models for degra- 
dation processes having more than one step, each 
with its own rate constant, it is, in general, nec- 
essary to have adequate models for the important 
individual steps. For example, when the individual 
processes can be observed, it may be possible to es- 
timate the effect that temperature (or other accel- 
erating variables) has on each of the rate constants. 

5.6 Other Units for Activation Energy 

The discussion and examples of the Arrhenius and 
Eyring relationships in Sections 5.2-5.4 used units of 
electronvolt for Ea and electronvolt per kelvin for k. 
These units for the Arrhenius model are used most 
commonly in applications involving electronics. In 



other areas of application (e.g., degradation of or- 
ganic materials such as paints and coatings, plastics, 
food and pharmaceuticals), it is more common to see 
Boltzmann's constant k in units of electronvolt re- 
placed with the universal gas constant in other units. 
For example, the gas constant is commonly given in 
units of kilojoule per mole kelvin [i.e., R = 8.31447 
kJ/(mol-K)]. In this case, Ea is activation energy 
in units of kilojoule per mole (kJ/mol). The corre- 
sponding Arrhenius acceleration factor is 

^jr(temp, tempf/, Ea) 

\^ ( 120.27 120.27 \1 

= exp Ea 

[ Vtemp^K tempK/_ 

The universal gas constant can also be expressed 
in units of kilocalorie per mole kelvin, kcal/ (mol-K) 
[i.e., R = 1.98588 kcal/(mol-K)]. In this case, Ea is 
in units of kilocalorie per mole (kcal/mol). The cor- 
responding Arrhenius acceleration factor is 

^jr(temp, tempj/, Ea) 

\^ ( 503.56 503.56 \1 

= exp Ea - . 

I VtempfjK tempK/_ 

It is also possible to use units of kJ/(mol-K) and 
kcal/(mol-K) for the Ea coefficient in the Eyring 
model. 

Although k is standard notation for Boltzmann's 
constant and R is standard notation for the univer- 
sal gas constant, we use k to denote either of these 
in the Arrhenius relationship. 

5.7 Temperature Cycling 

Some failure modes are caused by temperature 
cycling. In particular, temperature cycling causes 
thermal expansion and contraction which can induce 
mechanical stresses. Some failure modes caused by 
thermal cycling include: 

• Power on/off cycling of electronic equipment can 
damage integrated circuit encapsulement and sol- 
der joints. 

• Heat generated by take-off power-thrust in jet en- 
gines can cause crack initiation and growth in fan 
disks. 

• Power-up/power-down cycles can cause cracks in 
nuclear power plant heat exchanger tubes and tur- 
bine generator components. 

• Temperature cycling can lead to delamination in 
Inkjet printhead components. 
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As in fatigue testing, it is possible to accelerate 
thermal cycling failure modes by increasing either 
the frequency or amplitude of the cycles (increasing 
amplitude generally increases mechanical stress) . The 
most commonly used model for acceleration of ther- 
mal cycling is the Coffin-Manson relationship which 
says that the number of cycles to failure is 



(Atemp)ft 



where Atemp is the temperature range and 6 and 
Pi are properties of the material and test setup. 
This power-rule relationship explains the effect that 
temperature range has on the thermal-fatigue life 
cycles-to-failure distribution. Nelson (1990, page 86) 
suggests that for some metals, /3i ~ 2 and that for 
plastic encapsulements used for integrated circuits, 
/?! ~ 5. The Cofhn-Manson relationship was origi- 
nally developed as an empirical model to describe 
the effect of temperature cycling on the failure of 
components in the hot part of a jet engine. See Nel- 
son (1990, page 86) for further discussion and refer- 
ences. 

Letting T be the random number of cycles to fail- 
ure (e.g., T = Ne where e is a random variable). The 
acceleration factor at Atemp, relative to Atemp^, is 



^J^(Atemp) 



T(Atemp^) 
T(Atemp) 



/ Atemp 
V Atempf; / 



There may be a Atemp threshold below which little 
or no fatigue damage is done during thermal cycling. 

Empirical evidence has shown that the effect of 
temperature cycling can depend importantly on 
tempj^g^^K, the maximum temperature in the cy- 
cling (e.g., if tempjjjg^^K is more than 0.2 or 0.3 times 
a metal's melting point). The cycles-to-failure distri- 
bution for temperature cycling can also depend on 
the cycling rate (e.g., due to heat buildup). An em- 
pirical extension of the Coffin-Manson relationship 
that describes such dependencies is 



N: 



1 



(Atemp)^i (freq)^2 



X exp 



fEa X 11605 \ 



tempj, 



where freq is the cycling frequency and Ea is a 
quasi-activation energy. 

As with all acceleration models, caution must be 
used when using such a model outside the range of 
available data and past experience. 



6. USING HUMIDITY TO ACCELERATE 
REACTION RATES 

Humidity is another commonly used accelerating 
variable, particularly for failure mechanisms involv- 
ing corrosion and certain kinds of chemical degrada- 
tion. 

Example 7 (Accelerated life test of a printed wiring 
board). Figure 6 shows data from an ALT of printed 
circuit boards, illustrating the use of humidity as an 
accelerating variable. This is a subset of the larger 
experiment described by LuValle, Welsher and Mitchell 
(1986), involving acceleration with temperature, hu- 
midity and voltage. A table of the data is given in 
Meeker and LuValle (1995) and in Meeker and Es- 
cobar (1998). Figure 6 shows clearly that failures 
occur earlier at higher levels of humidity. 

Vapor density measures the amount of water va- 
por in a volume of air in units of mass per unit vol- 
ume. Partial vapor pressure (sometimes simply re- 
ferred to as "vapor pressure" ) is closely related and 
measures that part of the total air pressure exerted 
by the water molecules in the air. Partial vapor pres- 
sure is approximately proportional to vapor density. 
The partial vapor pressure at which molecules are 
evaporating and condensing from the surface of wa- 
ter at the same rate is the saturation vapor pressure. 
For a fixed amount of moisture in the air, saturation 
vapor pressure increases with temperature. 

Relative humidity is usually defined as 

Vapor Pressure 



RH: 
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Fig. 6. Scatterplot of printed circuit board accelerated life 
test data. Censored observations are indicated by A. There 
are 48 censored observations at 4078 hours in the 49.5% RH 
test and 11 censored observations at 3067 hours in the 62.8% 
RH test. 
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and is commonly expressed percent. For most 
failure mechanisms, physical/chemical theory sug- 
gests that RH is the appropriate scale in which to 
relate reaction rate to humidity especially if temper- 
ature is also to be used as an accelerating variable 
(Klinger, 1991). 

A variety of different humidity models (mostly 
empirical but a few with some physical basis) have 
been suggested for different kinds of failure mech- 
anisms. Much of this work has been motivated by 
concerns about the effect of environmental humidity 
on plastic-packaged electronic devices. Humidity is 
also an important factor in the service-life distribu- 
tion of paints and coatings. In most test applications 
where humidity is used as an accelerating variable, it 
is used in conjunction with temperature. For exam- 
ple. Peck (1986) presents data and models relating 
life of semiconductor electronic components to hu- 
midity and temperature. See also Peck and Zierdt 
(1974) and Joyce et al. (1985). Gillen and Mead 
(1980) describe a kinetic approach for modeling ac- 
celerated aging data. LuValle, Welsher and Mitchell 
(1986) describe the analysis of time-to-failure data 
on printed circuit boards that have been tested at 
higher than usual temperature, humidity and volt- 
age. They suggest ALT models based on the physics 
of failure. Chapter 2 of Nelson (1990) and Boccaletti 
et al. (1989) review and compare a number of dif- 
ferent humidity models. 

The Eyring/Arrhenius temperature-humidity ac- 
celeration relationship in the form of (14) uses xi = 
11605/temp K, X2 = log(RH) and X3 = X1X2 where 
RH is relative humidity, expressed as a proportion. 
An alternative humidity relationship suggested by 
Klinger (1991), on the basis of a simple kinetic model 
for corrosion, uses the term X2 = log[RH/(l — RH)] (a 
logistic transformation) instead. 

In most applications where it is used as an accel- 
erating variable, higher humidity increases degrada- 
tion rates and leads to earlier failures. In applica- 
tions where drying is the failure mechanism, how- 
ever, an artificial environment with lower humidity 
can be used to accelerate a test. 

7. ACCELERATION MODEL FOR 
PHOTODEGRADATION 

Many organic compounds degrade chemically when 
exposed to ultraviolet (UV) radiation. Such degra- 
dation is known as photodegradation. This section 
describes models that have been used to study pho- 
todegradation and that are useful when analyzing 



data from accelerated photodegradation tests. Many 
of the ideas in this section originated from early re- 
search into the effects of light on photographic emul- 
sions (e.g., James, 1977) and the effect that UV ex- 
posure has on causing skin cancer (e.g., Blum, 1959). 
Important applications include prediction of service 
life of products exposed to UV radiation (outdoor 
weathering) and fiber-optic systems. 

7.1 Time Scale and Model for Total Effective 
UV Dosage 

As described in Martin et al. (1996), the appro- 
priate time scale for photodegradation is the total 
(i.e., cumulative) effective UV dosage, denoted by 
DTot- Intuitively, this total effective dosage can be 
thought of as the cumulative number of photons ab- 
sorbed into the degrading material and that cause 
chemical change. The total effective UV dosage at 
real time t can be expressed as 

(8) DTot(t)= /*Dinst (r)dT 

Jo 

where the instantaneous effective UV dosage at real 
time T is 

Dinst (t) = / Dinst (t, A) dX 

(9) =(^'e^{\t) 

X {l-exp[-A{X)]}<P{X)dX. 

Here Eo(X, r) is the spectral irradiance (or intensity) 
of the light source at time r (both artificial and nat- 
ural light sources have potentially time-dependent 
mixtures of light at different wavelengths, denoted 
by A), [1 — exp(— A(A))] is the spectral absorbance of 
the material being exposed (damage is caused only 
by photons that are absorbed into the material) , and 
(f){X) is a quasi-quantum efficiency of the absorbed 
radiation (allowing for the fact that photons at cer- 
tain wavelengths have a higher probability of caus- 
ing damage than others) . The functions £"0 and A in 
the integrand of (9) can either be measured directly 
or estimated from data and the function (j){X) can 
be estimated from data. A simple log-linear model is 
commonly used to describe quasi-quantum efficiency 
as a function of wavelength. That is, 

</.(A)=exp(/3o + /3iA). 

The integrals over wavelength, like that in (9), are 
typically taken over the UV-B band (290 nm to 
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320 nm), as this is the range of wavelengths over 
which both (p{X) and Eo(X,t) are importantly dif- 
ferent from 0. Longer wavelengths (in the UV-A 
band) are not terribly harmful to organic materials 
[(/)(A) ^ 0] . Shorter wavelengths (in the UV-C band) 
have more energy, but are generally filtered out by 
ozone in the atmosphere [EQ{X,t) ~0]. 

7.2 Additivity 

Implicit in the model (9) is the assumption of ad- 
ditivity. Additivity implies, in this setting, that the 
photo-effectiveness of a source is equal to the sum 
of the effectiveness of its spectral components. This 
part of the model makes it relatively easy to con- 
duct exposure tests with specific combinations of 
wavelengths [e.g., by using selected band-pass fil- 
ters to define Eq{X,t) functions as levels of spectral 
intensity in an experiment] to estimate the quasi- 
quantum efficiency function of A. Then the to- 
tal dosage model in (9) can be used to predict pho- 
todegradation under other combinations of wave- 
lengths [i.e., for other Eq{X,t) functions]. 

7.3 Reciprocity and Reciprocity Breakdown 

The intuitive idea behind reciprocity in photodegra- 
dation is that the time to reach a certain level of 
degradation is inversely proportional to the rate at 
which photons attack the material being degraded. 
Reciprocity breakdown occurs when the coefficient 
of proportionality changes with light intensity. Al- 
though reciprocity provides an adequate model for 
some degradation processes (particularly when the 
range of intensities used in experimentation and ac- 
tual applications is not too broad), some examples 
have been reported in which there is reciprocity break- 
down (e.g., Blum, 1959, and James, 1977). 

Light intensity can be affected by filters. Sunlight 
is filtered by the earth's atmosphere. In laboratory 
experiments, different neutral density filters are used 
to reduce the amount of light passing to specimens 
(without having an important effect on the wave- 
length spectrum), providing an assessment of the 
degree of reciprocity breakdown. Reciprocity implies 
that the effective time of exposure is 



a!(t) = CFxDTot(t) 
CF X 



JAi 



Dinst(T, A)dA(ir 



to concentrate light to achieve, say, "5 Suns" ac- 
celeration or CF = 5. A 50% neutral density filter 
in a laboratory experiment will provide deceleration 
corresponding to CF = 0.5. 

When there is evidence of reciprocity breakdown, 
the effective time of exposure is often modeled, em- 
pirically, by 



(10) 



d(t) = (CF)PxDTot(t) 



(CF)^ 



t l-\2 



J\ 



Dinst(T, A) dX dr 



where CF is an "acceleration factor." For example, 
commercial outdoor test exposure sites use mirrors 



Model (10) has been shown to fit data well and ex- 
perimental work in the photographic literature sug- 
gests that when there is reciprocity breakdown, the 
value of p does not depend strongly, if at all, on the 
wavelength A. A statistical test for p = 1 can be used 
to assess the reciprocity assumption. 

7.4 Model for Photodegradation and UV 
Intensity 

Degradation (or damage) D{t) at time t depends 
on environmental variables like UV, temp and RH, 
that may vary over time, say according to a mul- 
tivariable profile ^(t) = [UV,temp,RH, . . .]. Labora- 
tory tests are conducted in well-controlled environ- 
ments, usually holding these variables constant (al- 
though sometimes such variables are purposely changed 
during an experiment, as in step-stress accelerated 
tests). Interest often centers, however, on life in a 
variable environment. Figure 7 shows some typical 
sample paths (for FTIR peak at 1510 cm~^, repre- 
senting benzene ring mass loss) for several specimens 
of an epoxy exposed to UV radiation using a band- 
pass filter with a nominal center at 306 nm. Separate 
paths are shown for each combination of (10, 40, 60, 
100)% neutral density filters and 45 °C and 55 °C, 
as a function of total (cumulative) absorbed UV-B 
dosage. These sample paths might be modeled by a 
given functional form, 

D{t)=g{z), z = log[d(t)]-//, 

where z is scaled time and g{z) would usually be 
suggested by knowledge of the kinetic model (e.g., 
linear for zeroth-order kinetics and exponential for 
first-order kinetics), although empirical curve fitting 
may be adequate for purposes where the amount of 
extrapolation in the time dimension is not large. As 
in SAFT models, can be modeled as a function 
of explanatory variables like temperature and hu- 
midity when these variables affect the degradation 
rate. 
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8. VOLTAGE AND VOLTAGE-STRESS 
ACCELERATION 

Increasing voltage or voltage stress (electric field) 
is another commonly used method to accelerate fail- 
ure of electrical materials and components like light 
bulbs, capacitors, transformers, heaters and insula- 
tion. 

Voltage quantifies the amount of force needed to 
move an electric charge between two points. Phys- 
ically, voltage can be thought of as the amount of 
pressure behind an electrical current. Voltage stress 
quantifies voltage per unit of thickness across a di- 
electric and is measured in units of volt /thickness 
(e.g., V/mm or kV/mm). 

8.1 Voltage Acceleration Mechanisms 

Depending on the failure mode, higher voltage 
stress can: 

• accelerate failure-causing electrochemical reactions 
or the growth of failure-causing discontinuities in 
the dielectric material. 

• increase the voltage stress relative to dielectric 
strength of a specimen. Units at higher stress will 
tend to fail sooner than those at lower stress. 

Sometimes one or the other of these effects will be 
the primary cause of failure. In other cases, both 
effects will be important. 

Example 8 (Accelerated life test of insulation for 
generator armature bars). Doganaksoy, Hahn and 



Meeker (2003) discuss an ALT for a new mica-based 
insulation design for generator armature bars (GABs). 
Degradation of an organic binder in the insulation 
causes a decrease in voltage strength and this was 
the primary cause of failure in the insulation. The 
insulation was designed for use at a voltage stress of 
120 V/mm. Voltage-endurance tests were conducted 
on 15 electrodes at each of five accelerated voltage 
levels between 170 V/mm and 220 V/mm (i.e., a 
total of 75 electrodes). Each test was run for 6480 
hours at which point 39 of the electrodes had not 
yet failed. Table 1 gives the data from these tests. 
The insulation engineers were interested in the 0.01 
and 0.05 quantiles of lifetime at the use condition of 
120 V/mm. Figure 8 plots the insulation lifetimes 
against voltage stress. 

8.2 Inverse Power Relationship 

The inverse power relationship is frequently used 
to describe the effect that stresses like voltage and 
pressure have on lifetime. Voltage is used in the fol- 
lowing discussion. When the thickness of a dielectric 
material or insulation is constant, voltage is propor- 
tional to voltage stress. Let volt denote voltage and 
let volt;7 be the voltage at use conditions. The life- 
time at stress level volt is given by 



T(volt) 



r(volt 
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spectrum] . 
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8. GAB insulation data. Scatterplot of life versus volt- 
Censored observations are indicated by A. 



where in general, is negative. The model has 
SAFT form with acceleration factor 



^J^(volt) 



(11) 



^^(volt, voltij,l3i] 
r(voltc/) 



T(volt 
volt 



-/3i 



,VOlt[/^ 

If r(volt[/) has a log-location-scale distribution with 
parameters fiu and a, then T(volt) also has a log- 
location-scale distribution with fi = Pq + Pix, where 

XU = log(voltc/), X = log(volt), /3o = /if/ - Pixu 

and a does not depend on x. 

Example 9 {Time acceleration for GAB insu- 
lation). For the GAB insulation data in Exam- 
ple 8, an estimate for f3i is Pi = —9 (methods for 
computing such estimates are described in Meeker 
and Escobar, 1998, Chapter 19). Recall that the de- 
sign voltage stress is volt[/ = 120 V/mm and con- 
sider testing at volt = 170 V/mm. Thus, using Pi = 
Pi, ^J^(170) = (170/120)9 w 23. Thus by increasing 




MO 1SS 130 13S 1JI> 145 



T55 165 



Fig. 9. Time- acceleration factor as a function of voltage 
stress and exponent —f3i — —7, —9, —11. 



voltage stress from 120 V/mm to 170 V/mm, one 
estimates that lifetime is shortened by a factor of 
l/AT{170) « 1/23 = 0.04. Figure 9 plots AJ^ ver- 
sus volt for Pi = —7, —9, —11. Using direct compu- 
tations or from the plot, one obtains AJ-{170) ~ 11 
for Pi = -7 and ^J^(170) « 46 for Pi = -11. 

Example 10 (Accelerated life test of a mylar- 
polyurethane insulation). Meeker and Escobar (1998, 
Section 19.3) reanalyzed ALT data from a special 
type of mylar-polyurethane insulation used in high- 
performance electromagnets. The data, originally from 
Kalkanis and Rosso (1989), give time to dielectric 
breakdown of units tested at (100.3, 122.4, 157.1, 
219.0, 361.4) kV/mm. The purpose of the ALT was 
to evaluate the reliability of the insulating struc- 
ture and to estimate the life distribution at sys- 
tem design voltages, assumed to be 50 kV/mm. Fig- 
ure 10 shows that failures occur much sooner at high 
voltage stress. Except for the 361.4 kV/mm data, 
the relationship between log life and log voltage ap- 



Table 1 
GAB insulation data 



Voltage stress Lifetime 
(V/mm) (thousand hours) 



170 15 censored" 

190 3.248, 4.052, 5.304, 12 censored" 

200 1.759, 3.645, 3.706, 3.726, 3.990, 5.153, 6.368, 8 censored" 

210 1.401, 2.829, 2.941, 2.991, 3.311, 3.364, 3.474, 4.902, 5.639, 6.021, 6.456, 4 censored" 

220 0.401, 1.297, 1.342, 1.999, 2.075, 2.196, 2.885, 3.019, 3.550, 3.566, 3.610, 3.659, 3.687, 4.152, 5.572 



"Units were censored at 6.480 thousand hours. 
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Fig. 10. Inverse power relationship-lognormal model fit- 
ted to the mylar-polyurethane data (also showing the 
361.4 kV/mm data omitted from the ML estimation). 

pears to be approximately linear. Meeker and Es- 
cobar (1998), in their reanalysis, omitted the 361.4 
kV/mm data because it is clear that a new failure 
mode had manifested itself at this highest level of 
voltage stress. Insulation engineers have suggested 
to us that the new failure mode was likely caused 
by thermal buildup that was not important at lower 
levels of voltage stress. 

8.3 Physical Motivation for the Inverse Power 
Relationship for Voltage-Stress Acceleration 

The inverse power relationship is widely used to 
model life as a function of pressure-like accelerat- 
ing variables (e.g., stress, pressure, voltage stress). 
This relationship is generally considered to be an 
empirical model because it has no formal basis from 
knowledge of the physics/chemistry of the modeled 
failure modes. It is commonly used because engi- 
neers have found, over time, that it often provides a 
useful description of certain kinds of AT data. 

This section presents a simple physical motivation 
for the inverse power relationship for voltage-stress 
acceleration under constant temperature situations. 
Section 9.2 describes a more general model for volt- 
age acceleration involving a combination of temper- 
ature and voltage acceleration. 

This discussion here is for insulation. The ideas 
extend, however, to other dielectric materials, prod- 
ucts and devices like insulating fluids, transform- 
ers, capacitors, adhesives, conduits and containers 
that can be modeled by a stress-strength interfer- 
ence model. 

In applications, an insulator should not conduct 
an electrical current. An insulator has a character- 
istic dielectric strength which can be expected to be 
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Fig. 11. Dielectric strength degrading over time, relative to 
voltage-stress levels (horizontal lines). 

random from unit to unit. The dielectric strength 
of an insulation specimen operating in a specific en- 
vironment at a specific voltage may degrade with 
time. Figure 11 shows a family of simple curves to 
model degradation and unit-to-unit variability in di- 
electric strength over time. The unit-to-unit vari- 
ability could be caused, for example, by materials or 
manufacturing variability. The horizontal lines rep- 
resent voltage-stress levels that might be present in 
actual operation or in an accelerated test. When a 
specimen's dielectric strength falls below the applied 
voltage stress, there will be flash-over, a short cir- 
cuit, or other failure-causing damage to the insula- 
tion. Analytically, suppose that degrading dielectric 
strength at age t can be expressed as 

V{t) = 5Q X t^/^i. 

Here, as in Section 5.4, failure occurs when ^{t) 
crosses Pf, the applied voltage stress, denoted by 
volt. In Figure 11, the unit-to-unit variability is 
in the 5q parameter. Equating T>{T) to volt and 
solving for failure time T gives 



Then the acceleration factor for volt versus voltf/ 
is 

.AJ^(volt) = ^J'(volt, volte/, /3i) 

_ T(voltt/) 
r(volt) 

_ / volt \ 
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which is an inverse power relationship, as in (11). 

To extend this model, suppose that higher voltage 
also leads to an increase in the degradation rate and 
that this increase is described with the degradation 
model 



V{t) = 6o[n{volt) X t] 



1/71 



where 



7^(volt) = 7o exp[72 log(volt)]. 



Suppose failure occurs when crosses 2?f , the ap- 
plied voltage stress, denoted by volt. Then equating 
ViT) to volt and solving for failure time T gives 
the failure time 

^/ N 1 /volt 

r(volt) = — 

^ ' 7^(volt) V So 

Then the ratio of failure times at volt^/ versus volt 
is the acceleration factor 



71 



^j^(volt) 



r(voit 



u) 



r(voit) V'voit[/ 



volt 



72-71 



which is again an inverse power relationship with 

Pi =71 -72- 

This motivation for the inverse power relationship 
described here is not based on any fundamental un- 
derstanding of what happens to the insulating ma- 
terial at the molecular level over time. As we de- 
scribe in Section 11, the use of such fundamental 
understanding could provide a better, more credible 
model for extrapolation. 

8.4 Other Inverse Power Relationships 

The inverse power relationship is also commonly 
used for other accelerating variables including pres- 
sure, cycling rate, electric current, stress and humid- 
ity. Some examples are given in Section 9. 

8.5 A More General Empirical Power 
Relationship: Box-Cox Transformations 

As shown in Section 8.2, the inverse power rela- 
tionship induces a log-transformation in volt giving 
the model = /3o + Pix, where x = log(volt). There 
might be other transformations of volt that could 
provide a better description of the data. A general, 
and useful, approach is to expand the formulation of 
the model by adding a parameter or parameters and 
investigating the effect of perturbing the added pa- 
rameter (s), to see the effect on answers to questions 
of interest. Here this approach is used to expand the 
inverse power relationship model. 



Suppose that Xi is a positive accelerating variable 
and X2 is a collection of other explanatory variables, 
some of which might be accelerating variables. Con- 
sider the model fj, = (3o + PiXi + where the 
/9's are unknown parameters. We start by replacing 
Xi with the more general Box-Cox transformation 
(Box and Cox, 1964) on Xi. In particular, we fit the 
model 

^ = 70 + 71 ^1 + 72^2 
where the 7's are unknown parameters and 
X^ - 1 



(12) 



A 



log(Xi), A = 0. 



The Box-Cox transformation (Box and Cox, 1964) 
was originally proposed as a simplifying transforma- 
tion for a response variable. Transformation of accel- 
erating and explanatory variables, however, provides 
a convenient extension of the accelerating modeling 
choices. The Box-Cox transformation includes all 
the power transformations and because Wi is a con- 
tinuous function of A, (12) provides a continuum of 
transformations for possible evaluation and model 
assessment. The Box-Cox transformation parame- 
ter A can be varied over some range of values (e.g., 
— 1 to 2) to see the effect of different voltage-life 
relationships on the fitted model and inferences of 
interest. The results from the analysis can be dis- 
played in a number of different ways. 

For fixed X2, the Box-Cox transformation model 
acceleration factor is 



^.Fbc(^i) 



exp|^ 



71 



if A / 0, 
if A = 0, 



where Xm are use conditions for the Xi accelerating 
variable. AJ- Bci^i) is monotone increasing in Xi if 
7i < and monotone decreasing in Xi if 71 > 0. 

Example 11 {Spring life test data). Meeker, Es- 
cobar and Zayac (2003) analyze spring accelerated 
life test data. Time is in units of kilocycles to failure. 
The explanatory variables are processing tempera- 
ture (Temp) in degrees Fahrenheit, spring compres- 
sion displacement (Stroke) in mils, and the categor- 
ical variable Method which takes the values New or 
Old. Springs that had not failed after 5000 kilocy- 
cles were coded as "Suspended." At the condition 50 
mils, 500 °F and the New processing method, there 
were no failures before 5000 kilocycles. All of the 



ACCELERATED TEST MODELS 



19 



««• VMul SMw Twin ftw. I 




^1,0 -as 0-0 0.5 1.0 1.5 £.0 
Slrofcfl Ben-Cox Tfansfomnallort Power 

Fig. 12. Plot of the ML estimate of the 0.10 quantile of 
spring life at 20 mils, 600 ° F, using the new method versus 
the Stroke displacement Box-Cox transformation parameter A 
with 95% confidence limits. 

other conditions had at least some failures, and at 
five of the twelve conditions all of the springs failed. 
At some of the conditions, one or more of the springs 
had not failed after 5000 kilocycles. 

Figure 12 (see Meeker, Escobar and Zayac, 2003) 
is a plot of the 0.10 Weibull quantile estimates versus 
A from —1 and 2. Approximate confidence intervals 
are also given. The plot illustrates the sensitivity 
of the 0.10 quantile estimate to the Box-Cox trans- 
formation. Note that A = corresponds to the log- 
transformation that is commonly used in fatigue life 
versus stress models. Also, A = 1 corresponds to no 
transformation (or, more precisely, a linear transfor- 
mation that affects the regression parameter values 
but not the underlying structure of the model) . Fig- 
ure 12 shows that fatigue life decreases by more than 
an order of magnitude as A moves from to 1. In 
particular, the ML estimate of the 0.10 quantile de- 
creases from 900 megacycles to 84 megacycles when 
A is changed from to 1. 

Figure 13 is a profile likelihood plot for the Box- 
Cox A parameter, providing a visualization of what 
the data say about the value of this parameter. In 
this case the peak is at a value of A close to 0; 
this is in agreement with the commonly used fa- 
tigue life/stress model. Values of A close to 1 are 
less plausible, but cannot be ruled out, based on 
these data alone. The engineers, based on experi- 
ence with the same failure mode and similar mate- 
rials, felt that the actual value of A was near (cor- 
responding to the log-transformation) and almost 



certainly less than 1. Thus a conservative decision 
could be made by designing with an assumed value 
of A = 1 . Even the somewhat optimistic evaluation 
using A = would not meet the 500 megacycle tar- 
get life. 

Meeker, Escobar and Zayac (2003) also discuss the 
sensitivity to the assumed form of the temperature- 
life relationship and the sensitivity to changes in the 
assumed distribution. 

9. ACCELERATION MODELS WITH MORE 
THAN ONE ACCELERATING VARIABLE 

Some accelerated tests use more than one acceler- 
ating variable. Such tests might be suggested when 
it is known that two or more potential accelerat- 
ing variables contribute to degradation and failure. 
Using two or more variables may provide needed 
time-acceleration without requiring levels of the in- 
dividual accelerating variables to be too high. Some 
accelerated tests include engineering variables that 
are not accelerating variables. Examples include ma- 
terial type, design, operator, and so on. 

9.1 Generalized Eyring Relationship 

The generalized Eyring relationship extends the 
Eyring relationship in Section 5.3, allowing for one 
or more nonthermal accelerating variables (such as 
humidity or voltage). For one additional nonther- 
mal accelerating variable X, the model, in terms of 
reaction rate, can be written as 



7e(temp,X) 
(13) = 70 X (temp K)"" x exp 



-71 



k X temp K 
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Fig. 13. Profile UkeUhood plot for the Stroke Box-Cox trans- 
formation parameter A m the spring life model. 
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X exp 72X + 



73X 



k X temp K 



where X is a function of the nonthermal stress. The 
parameters 71 = Ea (activation energy) and 70, 72, 
73 are characteristics of the particular physical/ 
chemical process. Additional factors like the one on 
the right-hand side of (13) can be added for other 
nonthermal accelerating variables. 

In the following sections, following common prac- 
tice, we set (tempK)"^ = 1, using what is essentially 
the Arrhenius temperature-acceleration relationship. 
These sections describe some important special-case 
applications of this more general model. If the un- 
derlying model relating the degradation process to 
failure is a SAFT model, then, as in Section 5.2, 
the generalized Eyring relationship can be used to 
describe the relationship between times at different 
sets of conditions temp and X. In particular, the 
acceleration factor relative to use conditions temp^ 
and Xu is 



AJ='{temp, X) 



7^(temp,X) 
TZ{tem-p^,Xu)' 



The same approach used in Section 5.4 shows the 
effect of accelerating variables on time to failure. 
For example, suppose that T(temp^) (time at use or 
some other baseline temperature) has a log-location- 
scale distribution with parameters fiu and a. Then 
T(temp) has the same log- location-scale distribution 
with 



(14) 



fi = - log[^J^(temp, X)] 

= Po+ PlXi + (52X2 + /33X1X2 



data by including second-order terms for both accel- 
erating variables. 

To derive the time-acceleration factor for the ex- 
tended Arrhenius relationship with temp and volt, 
one can follow steps analogous to those outlined in 
Section 8.2. Using (13) with X = log(volt), one ob- 
tains 



where (3i = Ea, P2 = -72, Ps = -73, xi = 11605/ 
(tempK), X2= X and Po = Hu - Pixw - P2X2U - 

P3XIUX2U- 

9.2 Temperature-Voltage Acceleration Models 

Many different models have been used to describe 
the effect of the combination of temperature and 
voltage on acceleration. For instance. Meeker and 
Escobar (1998, Section 17.7) analyzed data from a 
study relating voltage and temperature to the fail- 
ure of glass capacitors. They modeled the location 
parameter of log-lifetime as a simple linear function 
of temp °C and volt. The generalized Eyring rela- 
tionship in Section 13 can also be used with X = 
log(volt), as done in Boyko and Gerlach (1989). 
Klinger (1991) modeled the Boyko and Gerlach (1989) 



7^(temp, volt 

= 70 X exp 

X exp 



-Ea 



k X temp K 



72 log (volt) + 



73 log (volt) 
k X temp K 



Again, failure occurs when the dielectric strength 
crosses the applied voltage stress, that is, ^{t) = 
volt. This occurs at time 



T(temp, volt) 



1 



7^(temp, volt) 
From this, one computes 
AJ- (temp, volt) 

T {tempjj, volt If) 



voltyi 



r(temp, volt) 

■ exp[Ea{xiu - xi)] X 
X {exp[j;i log(volt) 



volt y2-7i 

VOlt(/ / 

xlf/log(volt{/)]}'^^ 



where xiu = 11605/(temp^K) and xi = 11605/ 
(tempK). When 73 = 0, there is no interaction be- 
tween temperature and voltage. In this case, 
^jF(temp, volt) can be factored into two terms, one 
that involves temperature only and another term 
that involves voltage only. Thus, if there is no in- 
teraction, the contribution of temperature (voltage) 
to acceleration is the same at all levels of voltage 
(levels of temperature). 

9.3 Temperature-Current Density 
Acceleration Models 

d'Heurle and Ho (1978) and Ghate (1982) studied 
the effect of increased current density (A/cm^) on 
electromigration in microelectronic aluminum con- 
ductors. High current densities cause atoms to move 
more rapidly, eventually causing extrusion or voids 
that lead to component failure. ATs for electromi- 
gration often use increased current density and tem- 
perature to accelerate the test. An extended Arrhe- 
nius relationship could be appropriate for such data. 
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In particular, when T has a log-location-scale distri- 
bution, then (13) applies with x\ = 11605/temp K, 
X2 = log(current). The model with /?3 = (without 
interaction) is known as "Black's equation" (Black, 
1969). 

Example 12 [Light emitting diode {LED) reli- 
ability]. A degradation study on a light emitting 



diode (LED) device was conducted to study the ef- 
fect of current and temperature on light output over 
time and to predict life at use conditions. A unit 
was said to have failed if its light output was re- 
duced to 60% of its initial value. Two levels of cur- 
rent and six levels of temperature were used in the 
test. Figure 14 shows the LED light output data ver- 




Hours 



114Cslsiu5;3amA 



low 4M0 



130CelsiiJs;40inA 




1000 4000 



Fig. 14. Relative change m light output from 138 hours at different levels of temperature and current. Relative change is in 
the linear scale and hours is in the square-root scale which linearizes the response as a function of time. 




Fig. 15. LED device data. Lognormal multiple probability plot showing the fitted Arrhenius-inverse power relationship log- 
normal model (with no interaction) for the failure LED data. The plots also show the estimate of F{t) at use conditions, 20 ° C 
and 20 mA. 
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sus time in hours, in the square-root scale. No units 
had failed during the test. For a simple method of 
degradation analysis, predicted pseudo failure times 
are obtained by using ordinary least squares to fit 
a line through each sample path on the square-root 
scale, for "Hours," and the linear scale for "Rela- 
tive Change." Figure 15 shows the ML fit of the 
Arrhenius- inverse power relationship lognormal model 
(with no interaction) for the pseudo failure LED 
data. The data at 130 °C and 40 mA were omit- 
ted in the model fitting because it was determined 
that a new failure mode had manifested itself at that 
highest level of test conditions (initial efforts by en- 
gineers to use the bad data had resulted in physi- 
cally impossible estimates of life at the use condi- 
tions). Figure 15 also shows the estimate of F{t) at 
use conditions of 20 °C and 20 mA. 

9.4 Temperature-Humidity Acceleration Models 

Relative humidity is another environmental vari- 
able that can be combined with temperature to ac- 
celerate corrosion or other chemical reactions. 
Examples of applications include paints and coat- 
ings, electronic devices and electronic semiconduc- 
tor parts, circuit boards, permalloy specimens, foods 
and pharmaceuticals. Although most ALT models 
that include humidity were derived empirically, some 
humidity models have a physical basis. For example, 
Gillen and Mead (1980) and Khnger (1991) stud- 
ied kinetic models relating aging with humidity. Lu- 
Valle, Welsher and Mitchell (1986) provided physical 
basis for studying the effect of humidity, tempera- 
ture and voltage on the failure of circuit boards. See 
Boccaletti et al. (1989), Chapter 2 of Nelson (1990), 
Joyce et al. (1985), Peck (1986) and Peck and Zierdt 
(1974) for ALT applications involving temperature 
and humidity. 

The extended Arrhenius relationship (13) applied 
to ALTs with temperature and humidity uses xi = 
11605/tenip K, X2 = log(RH) and X3 = X1X2 where 
RH is a proportion denoting relative humidity. The 
case when 03 = (no temperature-humidity interac- 
tion) is known as "Peck's relationship" and was used 
by Peck (1986) to study failures of epoxy packing. 
Klinger (1991) suggested the term X2 = log[RH/(l — 
RH)] instead of log(RH). This alternative relationship 
is based on a kinetic model for corrosion. 

9.5 Modeling Photodegradation Temperature 
and Humidity Effects 

When modeling photodegradation, as described in 
Section 7, it is often necessary to account for the ef- 
fect of temperature and humidity. The Arrhenius 



rate reaction model (3) can be used to scale time 
(or dosage) in the usual manner. Humidity is also 
known to affect photodegradation rate. Sometimes 
the rate of degradation will be directly affected by 
moisture content of the degrading material. In this 
case one can use a model such as described in Burch, 
Martin and VanLandingham (2002) to predict mois- 
ture content as a function of relative humidity. 

Combining these model terms with the log of total 
effective UV dosage from (8) gives 

log(d;CF,p) = log[DTot(t)] + P X log(CF), 

En 



k X tempK 



+ C X MC(RH), 



where tempK is temperature in kelvin, MC(RH) is a 
model prediction of moisture content, as a function 
of relative humidity, k is Boltzmann's constant, Ea 
is a quasi-activation energy, and (3q and C are pa- 
rameters that are characteristic of the material and 
the degradation process. Figure 16 shows the same 
data displayed in Figure 7, except that the time scale 
for the data has been adjusted for the differences in 
both the neutral density filters and the two different 
levels of temperature, bringing all data to the scale 
of a 100% neutral density filter and 55 °C. 

9.6 The Danger of Induced Interactions 

As illustrated in Example 13 of Pascual, Meeker 
and Escobar (2006) (using data related to Exam- 
ple 12 of this paper), interactions can cause difficulty 
in the interpretation of accelerated test results. Ini- 
tially, in that example, an interaction term had been 
used in the model fitting to provide a model that fits 
the data better at the high levels of temperature and 
current density that had been used in the test. Ex- 
trapolation to the use conditions, however, produced 
estimates of life that were shorter than the test con- 
ditions! The problem was that with the interaction 
term, there was a saddle point in the response sur- 
face, outside of the range of the data. Extrapolation 
beyond the saddle point resulted in nonsensical pre- 
dictions that lower temperature and current would 
lead to shorter life (in effect, the extrapolation was 
using a quadratic model). 

It is important to choose the definition of acceler- 
ated test experimental factors with care. Inappropri- 
ate choices can induce strong interactions between 
the factors. Suppose, for example, that there is no 
interaction between the factors voltage stress and 
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specimen thickness in an acceleration model for di- 
electric failure. Then 

^ = /3o + /?! Thickness + /32 Voltage Stress 

where thickness is measured in mm and Voltage 
Stress = Voltage/Thickness is measured in V/mm. 
If the model is written in terms of thickness and 
voltage, 

fJ- = /3o + /?! Thickness 

+ /32 Voltage Stress x Thickness. 

Thus, the variables voltage and thickness would have 
a strong interaction. 

Similarly, if RH and temperature have no interac- 
tion in the acceleration model for a corrosion failure 
mechanism, it is easy to show that the strong effect 
that temperature has on saturation vapor pressure 
would imply that the factors temperature and vapor 
pressure would have a strong interaction. 

These concerns are related to the "sliding factor" 
ideas described in Phadke (1989, Section 6.4). 

10. COMMENTS ON THE APPLICATION OF 
ACCELERATION MODELS 

10.1 Concerns About Extrapolation 

All accelerated tests involve extrapolation. The 
use of extrapolation in applications requires justi- 
fication. It is always best when the needed justifica- 
tion comes from detailed physical/chemical knowl- 
edge of the effect of the accelerating variable on the 



failure mechanism. As a practical matter, however, 
such knowledge is often lacking, as are time and re- 
sources for acquiring the needed knowledge. Empir- 
ical relationships are often used as justification, but 
rarely are data available to check the relationship 
over the entire range of interest for the accelerating 
variable(s). 

Evans (1977) makes the important point that the 
need to make rapid reliability assessments and the 
fact that accelerated tests may be "the only game 
in town" are not sufficient to justify the use of the 
method. Justification must be based on physical mod- 
els or empirical evidence. Evans (1991) describes dif- 
ficulties with accelerated testing and suggests the 
use of sensitivity analysis, such as that described 
in Meeker, Escobar and Zayac (2003). He also com- 
ments that acceleration factors of 10 "are not un- 
reasonable" but that "factors much larger than that 
tend to be figments of the imagination and lots of 
correct but irrelevant arithmetic." 

10.2 Some Basic Guidelines 

Some guidelines for the use of acceleration models 
include: 

• ATs must generate the same failure mode occur- 
ring in the field. 

Generally, accelerated tests are used to obtain 
information about one particular, relatively sim- 
ple failure mechanism (or corresponding degra- 
dation measure). If there is more than one fail- 
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ure mode, it is possible that the different failure 
mechanisms will be accelerated at different rates. 
Then, unless this is accounted for in the modeling 
and analysis, estimates could be seriously incor- 
rect when extrapolating to lower use levels of the 
accelerating variables. 

Accelerating variables should be chosen to corre- 
spond with variables that cause actual failures. 
It is useful to investigate previous attempts to 
accelerate failure mechanisms similar to the ones 
of interest. There are many research reports and 
papers that have been published in the physics 
of failure literature. The annual Proceedings of 
the International Reliability Physics Symposium, 
sponsored by the IEEE Electron Devices Society 
and the IEEE Reliability Society, contain numer- 
ous articles describing physical models for accel- 
eration and failure. 

Accelerated tests should be designed, as much as 
possible, to minimize the amount of extrapola- 
tion required, as described in Chapters 20 and 21 
of Meeker and Escobar (1998). High levels of ac- 
celerating variables can cause extraneous failure 
modes that would never occur at use levels of the 
accelerating variables. If extraneous failures are 
not recognized and properly handled, they can 
lead to seriously incorrect conclusions. Also, the 
relationship may not be accurate enough over a 
wide range of acceleration. 

In practice, it is difficult or impractical to verify 
acceleration relationships over the entire range of 
interest. Of course, accelerated test data should 
be used to look for departures from the assumed 
acceleration model. It is important to recognize, 
however, that the available data will generally pro- 
vide very little power to detect anything but the 
most serious model inadequacies. Typically there 
is no useful diagnostic information about possible 
model inadequacies at accelerating variable levels 
close to use conditions. 

Simple models with the right shape have generally 
proven to be more useful than elaborate multipa- 
rameter models. 

Sensitivity analyses should be used to assess the 
effect of perturbing uncertain inputs (e.g., inputs 
related to model assumptions). 
Accelerated test programs should be planned and 
conducted by teams including individuals knowl- 
edgeable about the product and its use environ- 
ment, the physical/chemical/mechanical aspects 
of the failure mode, and the statistical aspects of 
the design and analysis of reliability experiments. 



11. FUTURE RESEARCH IN THE 
DEVELOPMENT OF ACCELERATED TEST 
MODELS 

Research in the development of accelerated test 
models is a multidisciplinary activity. Statisticians 
have an important role to play on the teams of sci- 
entists that develop and use accelerated test models. 
On such a team engineers and scientists are primar- 
ily responsible for: 

• Identifying and enumerating possible failure modes 
and, for new products, predicting all possible life- 
limiting failure modes. 

• Understanding the physical/chemical failure mech- 
anisms that lead to a product's failure modes and 
for identifying accelerating variables that can be 
used to accelerate the failure mechanism. 

• Suggesting deterministic physical/chemical math- 
ematical relationships between the rate of the fail- 
ure mechanisms and the accelerating variable(s). 
When such a relationship is not available, they 
may be able to provide guidance from standard 
practice or previous experience with similar prod- 
ucts and materials. 

Statisticians are important for: 

• Planning appropriate experiments. Accelerated test 
programs often start with simple experiments to 
understand the failure modes that can occur and 
the behavior of mechanisms that can cause fail- 
ures. Use of traditional methods for designed ex- 
periments is important for these tests. In addition, 
there may be special features of accelerated tests 
(e.g., censoring) that require special test plan- 
ning methods (see Chapter 6 of Nelson, 1990, and 
Chapter 20 of Meeker and Escobar, 1998, for dis- 
cussion of accelerated test planning). 

• Providing expertise in the analysis of data aris- 
ing from preliminary studies and the accelerated 
tests themselves. Features such as censoring, mul- 
tiple failure modes and models that are nonlin- 
ear in the parameters are common. Methods for 
detecting model departures are particularly im- 
portant. Model departures may suggest problems 
with the data, sources of experimental variability 
that might be eliminated or problems with the 
suggested model that may suggest changes to the 
proposed model. 

• Identifying sources of variability in experimental 
data (either degradation data or life data) that re- 
flect actual variability in the failure mechanism. 
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For this purpose it is generally useful to compare 
field data. In most applications there will be ad- 
ditional variability in field data and it is impor- 
tant to understand the differences in order to de- 
sign appropriate laboratory experiments and to 
be able to draw useful conclusions and predictions 
from laboratory data. 
• Working within cross-disciplinary teams to de- 
velop statistical models of failure and accelera- 
tion based on fundamental understanding of fail- 
ure mechanisms. With fundamental understand- 
ing of failure mechanisms and knowledge of sources 
of variability, it is possible to develop, from first 
principles, failure-time distributions and acceler- 
ation models. Better, more credible, models for 
extrapolation should result from such modeling 
efforts if the assumptions and other inputs are ac- 
curate. Some examples of where such models have 
been developed include 

- Fisher and Tippett (1928) derived the asymp- 
totic distributions of extreme values and these 
results provide fundamental motivation for use 
of distributions such as the Weibull distribution 
(one of the three distributions of minima) in re- 
liability applications in which failure is caused 
by the first of many similar possible events (e.g., 
the failure of a large system with many similar 
possible failure points, none of which dominates 
in the failure process). 

- Tweedie (1956), in modeling electrophoretic mea- 
surements, used the distribution of first passage 
time of a Brownian motion with drift to derive 
the inverse Gaussian distribution. 

- Birnbaum and Saunders (1969), in modeling 
time to fracture from a fatigue crack growth 
process, derived a distribution that is today 
known as the Birnbaum-Saunders distribution. 
This distribution can be thought of as a discrete- 
time analog of the inverse Gaussian distribu- 
tion. 

- Meeker and LuValle (1995), using a kinetic model 
for the growth of conducting filaments, devel- 
oped a probability distribution, dependent on 
the level of relative humidity, to predict the fail- 
ure time of printed circuit boards. 

- Meeker, Escobar and Lu (1998), using a kinetic 
model for chemical degradation inside an elec- 
tronic amplifying device, developed a probabil- 
ity distribution which, when combined with the 
Arrhenius model, could be used to predict the 
failure time of the devices. 



- LuValle et al. (1998) and LuValle, Lefevre and 
Kannan (2004, pages 200-206) used physics- 
based models to describe degradation processes. 
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